The Transition in the Two-Dimensional Step Model: A Kosterlitz-Thouless Transition 
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Evidence for a Kosterlitz-Thouless transition in the 2D step model is obtained from Monte Carlo 
determinations of the helicity modulus. It is argued that the free energy of a single vortex at the 
center of the system depends logarithmically on the system size in spite of the fact that the spin 
interaction is not harmonic for small differences in the spin angles. We conclude that this is the 
reason for the Kosterlitz-Thouless transition in the 2D step model and that the harmonic spin 
interaction not is a necessary requirement. 
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The phase transition in two-dimensional (2D) XY 
models is known to take place through the vortex un- 
binding mechanism due to Kosterlitz and Thouless (KT) 
|jj . From the principles of universality one expects this 
transition to remain the same independent of details of 
the system as e.g. the underlying lattice structure. The 
precise spin interaction potential U(cj>), where <fi is the an- 
gle difference between neighboring spins, is not supposed 
to be essential either. U (<f>) is however required to be pe- 
riodic in 2tt and it seems always to have been presumed 
that the interaction in addition has to be harmonic for 
small 4>. The harmonicity for small (f> has to do with the 
energetics for vortex formation. With a harmonic poten- 
tial the energy for a single vortex in a L x L lattice goes 
as InL, and in the classical argument by Kosterlitz and 
Thouless this property is crucial for the transition. 

The subject of the present Letter is the 2D step model 
which is an XY model with a spin interaction that has no 
harmonic component. U ((f)) is instead a step-like function 



U(<p) — — J sign (cos (j>) 



(1) 



Since this potential is flat around <f> = there is no In in- 
dependence of the energy for a single vortex. This en- 
ergy is instead independent of system size and one would 
therefore expect a non-vanishing density of free vortices 
at all finite temperatures, and consequently no phase 
transition. Against that background the evidence from 
simulations for a transition were very intriguing. The 
first clear evidence for a transition was obtained from a 
Monte Carlo (MC) study of the susceptibility and the 
specific heat g]. The increase of the susceptibility with 
lattice size was considered to suggest a phase transition 
at T « 1.1. Later simulations also provided evidence that 
the transition actually is in the same universality class as 
the harmonic XY models |J. The similarity of the be- 
havior close to the transition of an harmonic XY model 
and the step model also led these authors to question 
the vortex unbinding as the mechanism behind the KT 
transition j^] . How could a transition driven by vortices 
be altogether insensitive to the very different energy cost 
for vortices in the two models? 



In the present Letter we address the question of the 
necessity of an harmonic spin potential for the KT tran- 
sition by examining the behavior of the 2D step model. 
We first calculate the helicity modulus and show that the 
behavior of this quantity gives strong support for a KT 
transition. We then demonstrate that the cost in free 
energy for a single vortex at the center of the system in 
fact goes as InL. It is this feature that stabilizes the low 
temperature phase. Finally, we refine the arguments to 
obtain quantitatively satisfactory estimates. 

The helicity modulus T, is a convenient quantity for 
the study of KT transitions due to its universal value 
2T/n at the transition and the known form of the 

approach to this universal value with L ||. The usual 
procedure in MC simulations is to determine T from a 
correlation function which involves some derivatives of 
U(4>) 0. Clearly, with a step-like potential the deriva- 
tives of the potential cannot be calculated and this ex- 
pression cannot be used. A way out is to instead start 
from the defining expression for the helicity modulus 

d 2 F 



T = 



OA 2 



(2) 



A=0 



and perform the simulations with fluctuating twist 
boundary conditions || . In these simulations one collects 
a histogram of the total twist -P(A). Since the probabil- 
ity for a certain twist is related to the free energy through 
P(A) cx e~ F( - A ^ T Eq. (§ becomes 

d 2 InP 



T = —T 



dA 2 



(3) 



A=0 



The simulations are done with twist variables in the 
two directions, A x and A y , which beside the spin vari- 
ables 0i are updated with the Metropolis algorithm. 
With Tij a unit vector between nearest neighbors and 
A = (A x , A y ), the Hamiltonian may be written 



H 



(ij) 



U[9t- 0* - -r l7 ■ A I = 



Ufa. 



To get a good acceptance ratio for the twist update 
it is necessary to make use of L different twist variables 
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in each direction, with A x — J2k=i , (and similarly 
in the y direction) where k is the column (row) num- 
ber. In our simulations, which for convenience were for a 
0(256) model, we used the potential £/(</>) = — J for 129 
angle differences [— 7r/2, 7t/2] and +J for the remaining 
127 p. This choice is not expected to be important for 
the transition properties, but gives a slight shift of the 
transition temperature as compared to the potential of 
Eq. (|l|). The length of the runs were typically 5 x 10 8 /L 
sweeps through the lattice. 

In Fig. [I] we show the histogram P(A) from Monte 
Carlo simulations at T = 0.05J. Since the histogram 
is peaked around A = the figure immediately gives 
evidence for a low temperature phase with a finite stiff- 
ness. To determine T we fit a quadratic curve to InP 
for |A/tt| < 1/3 and obtain T/T = 0.789(4) (where the 
given error is one standard deviation). Note that this 
is slightly larger than the universal value 2/ir « 0.637, 
which is required for a stable low temperature phase. 
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FIG. 1. Twist histogram obtained through MC simulations 
at T = 0.05. The solid line is obtained by fitting the data to 
a quadratic form which gives T/T — 0.789. 

An important feature of the step model is the gap in 
excitation energies; there are no excitations with energy 
< 2 J. At T -c J the system is therefore at all times in 
one of its numerous ground states which means that the 
histogram P(A) is independent of temperature. From 
Eq. (||) then follows a linear temperature dependence for 
T. This is in contrast to harmonic XY models for which 
T/J approaches unity in the low-temperature limit. 

The temperature dependence of T is shown in Fig. |^ 
for several system sizes together with the dashed line for 
the universal jump condition 2T/tt. We note that the 
curves start out linearly at low temperatures, become 
size-dependent at T/J ss 0.75, cross the universal line 
and then drop down to zero. Beside the unusual linear 
temperature-dependence at low T this behavior is just 
as in an ordinary harmonic XY model and therefore pre- 
cisely what one would expect for a KT transition. 

To determine the KT temperature we make use of the 
finite size dependence of T |6) . We follow the procedure 
in Ref. |L0| of first fitting our MC data for T from a nar- 
row temperature interval to second order polynomials in 
T, one for each L, and then fit the data to the expression 
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FIG. 2. The helicity modulus versus temperature for L — 8 
through 256. The rapid decrease down to zero directly follows 
the crossing of the data with the universal line, 2T/n (dashed 
line). This behavior is suggestive of a KT transition. Note 
that even though T — » at low temperatures, the system 
remains in the low temperature phase since T/T > 2/7T. 
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A + 2\nL 



(4) 



which amounts to adjusting Tkt and A to get the best 
possible fit. Using data for L > 16 we obtained Tkt = 
0.765(6). Fig. |H illustrates the good fit of the data at 
the transition temperature to the line from Eq. (||). We 
consider the above MC data to be strong evidence for a 
KT transition. This is in agreement with the conclusion 
in Ref. || that the step model is in the same universality 
class as the harmonic XY model. 
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FIG. 3. Determination of the transition temperature. The 
solid line is Eq. (Q) and the data points are the helicity mod- 
ulus from our simulations. 

We now propose an analytical analysis in order to un- 
derstand this behavior. We focus on the properties in the 
low temperature phase where the angular differences are 
restricted to the low energy region, \<fi\ < tt/2. A central 
idea in the present letter is to note that, while in the har- 
monic XY model the spinwave-vortex interaction only is 
a smooth perturbation, in the step model this interaction 
leads to a dramatic and crucial effect: while the energy 
of a system with a single vortex fixed in the center of the 
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system is finite, the free energy of this system grows as 
InP This is due to the change in entropy of spinwave 
fluctuations for the configuration with the fixed vortex, 
as compared to the vortex free case. Note that this is the 
entropy associated with a fixed vortex, not the positional 
entropy associated with a free vortex's variable location. 

To demonstrate the existence of this spinwave entropy 
we consider the configuration of spins in Fig. where 
we slightly reorganize the spins and delete the links in 
the radial direction. We will return below to the approx- 
imation involved in this step. 
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• • a) b) 

FIG. 4. Arrangement of links used to a) argue for the 
change in spinwave entropy for a vortex at a fixed position 
and b) estimate the helicity modulus. 

The condition for having a positive vortex in the center 
of the system in Fig. ^a is that the phase rotates by 2ir 
along each of the circles. We therefore introduce <& r = 

<f> along the circle with radius r, which is a sum of ss 
2irr values and let f2 r (<I> r ) denote the number of possible 
combinations of the </>:s at distance r as a function of 
<F r , which is defined only for <I> r = 2nn (with integer 
n). The probability for having a (positive) vortex is then 
determined by the product 

P vort«ll-^y- (5) 

r— 1 v ' 

The fraction O r (27r)/r2 r (0) may be calculated if we tem- 
porarily open up a closed path that makes up a circle. 
Since we want the same number of links along this path 
we need one more spin variable at one of the endpoints, 
which we take to be independent of the other endpoint. 
$ r then becomes a sum of « 2irr independent variables 
<\>. With the average of cj) being equal to zero and its 
variance given by a 2 the distribution of $> r for this open 
path becomes a Gaussian with width 2irra 2 : 

^(^)ocexp(-^-,) (6) 

We now make use of the fact that the number of possible 
configurations for the open and closed paths are the same 
if the spins at the endpoints of the open path are equal. 
Since this condition is equal to having $ r = 2im we 
conclude that VL r {2izn) = ft°P cn (27m). From Eqs. (§) 
and (H) the probability for a vortex then becomes 



Pvort = exp ( ~ 'J ~ eX P lnL ) ' 

for the entropy of a fixed vortex we obtain 5 V ort = 
lnP vort = — \ lni, and the free energy for a vortex at a 
fixed position in a system of size L finally becomes 

7T 

Pvort = -TSyort = T -~ In L. (7) 

For the harmonic XY model the well-known argument 
for the phase transition gives the free energy for having 
a vortex at any of the L 2 positions as 

AF = (ttJ -2T)lnL, 

and the transition takes place at T/J — 2 /it. In the step 
model the corresponding expression becomes 

AP = Pvort - 2Tlni = T - 2) InP 

where the temperature, at first sight only appears as a 
prefactor. However, there is a hidden temperature depen- 
dence in the variance a 2 . At low enough temperatures 
the <f>ij are restricted to the interval [— 7r/2, 7t/2] but with 
increasing temperature the <pij will more often take val- 
ues outside this interval, and a 2 will increase. Therefore, 
if AP is positive at low temperatures it will turn neg- 
ative at some finite temperature and this will give the 
transition. But if a 2 > ir/2 already at low temperatures 
there will be no transition. 

To see how the local restrictions on the angle differ- 
ences give rise to the non-zero helicity modulus we now 
turn to a rectangular geometry and consider the deter- 
mination of T from the distribution of and A y . One 
point with examining the distribution of the twist is to 
give predictions that are easy to compare with MC simu- 
lations. In the simplest approximation we again delete all 
links in the perpendicular direction, as in Fig. |Jb. For a 
single row the number of configurations consistent with a 
certain total twist becomes f2 row (A) cx exp (— A 2 /2Lc 2 ), 
in analogy with Eq. (|^) , and the number of configurations 
for the whole system with L rows becomes 

fi(A) = [n mw (A)] L cx exp (~A 2 /2a 2 ) . 

For the free energy we then arrive at P(A) = TA 2 /2a 2 
and with Eq. (g) the helicity modulus becomes 

T = T/a 2 . (8) 

In the absence of perpendicular links and at low tem- 
peratures, the (f>ij have a rectangular distribution, and 
from elementary integrals one finds a 2 = tt 2 /12 w 0.822. 
Through Eq. (§) this gives T/T w 1.22 which is about 
50% larger than T/T w 0.789 from Monte Carlo simula- 
tions, cf. Fig. |l|. A better estimate will be obtained below 
by including some links in the perpendicular direction. 

Comparing Eqs. (0) and (||) we find F vort (L) — T7rlni 
which is the same relation as in the harmonic model. This 
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shows that our two different calculations are equivalent 
which is a consequence of using the same approximation 
in both cases, i.e. neglecting all links perpendicular to 
the direction of interest. 

We now discuss the assumption used above, namely 
that the qualitatively correct behavior in a certain di- 
rection may be obtained even though one neglects the 
perpendicular links. This assumption is true only if the 
relative reduction of the number of allowed configura- 
tions obtained by introducing perpendicular links is es- 
sentially independent of the total twist in the direction 
of interest. We argue that this is a plausible assumption 
by considering two sets of configurations: 1) the set of 
all twist free configurations and 2) the set of configura- 
tions with a twist A x . There is then a transformation 
(fiij + A/L — ► cf>ij on all the horizontal links that trans- 
forms each member in the twist free set into a correspond- 
ing one in the twisted set. Since this transformation not 
affects the angle difference at the perpendicular links, it 
follows that the effect of the perpendicular links will be to 
exclude the same number of configurations in both these 
sets and this suggests that the relative reduction due to 
the perpendicular links will be independent of A. How- 
ever, this argument only serves to make our assumption 
a reasonable one; it is not conclusive. There is nothing 
that guarantees that the relative reduction of the allowed 
configurations (with \<j>ij \ < tt/2 for all horizontal links) 
will be the same for the two different sets. That a certain 
member of the twist free set is allowed doesn't imply that 
the corresponding member in the twisted set is allowed 
too. 



which together give an effective variance a 2 s ■ 
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FIG. 5. The configuration of links used to get quantita- 
tively reasonable values. The starting point is given in panel 
a) and after integrating out the upper and lower rows on each 
triple row one obtains the links shown in panel b). 

Even though we have argued that the perpendicular 
links may be neglected in a qualitative discussion, they 
have to be included in order to get reasonable quantita- 
tive estimates, since they do affect the variance a 2 . To 
include some perpendicular links we consider the config- 
uration in Fig. |^a, where the perpendicular links have 
been deleted at every third row only. The approach is 
then to integrate out the upper and lower rows of spins 
in each of these triple rows to give L/3 one-dimensional 
rows, cf. Fig. ||b. These integrations may be done analyt- 
ically or by Monte Carlo simulations on that geometry. 
The integrations also give correlations to the neighboring 
and next neighboring links, (</></>') and ((fxf)") (cf. Fig. @), 



*-i((l>)V 



+ (#') + (#"). (9) 



With L/3 rows, as in Fig. ||b, the expression for the he- 
licity modulus becomes T = T/(3a 2 s ). From our MC 
simulations on the geometry of Fig. ||a we get a 2 — 0.672, 
(<jxf>') = -0.120, and = -0.026. The first two of 

these numbers are easily obtained by integrating analyt- 
ically with symbolic software. With Eq. (^) this gives 
a 2 s = 0.526 and T/T = 0.633 which is less than 20% off 
T/T = 0.789 obtained above. 



FIG. 6. Illustration to the link-link correlations in Eq. (j9|). 

To conclude we have presented evidence from simula- 
tions that the 2D step model actually undergoes a KT 
transition. We have argued that the reason for the sta- 
bility of the low temperature phase against the formation 
of free vortices is the In L-dependence of the free energy 
for a vortex at a fixed position. From these results we are 
led to the conclusion that the harmonic spin interaction 
is not a necessary condition for a KT transition in a 2D 
spin model, and that the KT transition is a more general 
phenomenon than has so far been recognized. 
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